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Abstract 

It is shown here that in neutral ion gases the thermal energy transport 
can occur in the form of new types of thermal soliton waves. The solitons 
can form under a vanishing net heating function, and for a quadratic 
net heating. It is predicted that these solitons play an important role 
in a diversity of terrestrial and astrophysical phenomena. We claim that 
the reported soliton waves can be observed under ordinary laboratory 
conditions. 
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1 Introduction 

Compressive wave-front propagation is one of the major causes of energy trans- 
port in thermally conducting gases, and gives rise to a diverse range of physical 
phenomena (see [1-3] for further references). It is known that except for the 
case of subsonic wave motion, all types of associated temperature perturbation 
cause thermal instability in such systems. This behavior is known as the Field 
criterion [4] , which allows only highly subsonic wave motion in a thermally sta- 
ble medium. A direct consequence of the Field criterion is that the pressure 
variations remain small, hence can be neglected in the analysis of compressive 
thermal waves. The resulting waves formation is highly nonlinear, and depends 
on the heating function. Linearization and phase space analysis have been used 
to determine the nature of such waves, which have indicated that apart from 
travelling wave formations, there are steady wave fronts in such cases. 
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In this work we show that in thermally conducting ion gases, the compres- 
sive wave propagation can occur in the form of stable, shape preserving wave 
forms, known as soliton waves. Solitons can be identified as travelling wave 
solutions that retain their amplitude throughout their propagation, and have 
constant amplitude in the asymptotic limits. We show that these soliton waves 
are very general feature of thermally conducting gases, both in astrophysical 
and under normal laboratory conditions. These play an important role in the 
energy transport in the media. We state the conditions under which such soliton 
waves can form in a neutral ion gas, under a linear or a quadratic or a zero net 
heating function. 



2 Basic Equations 

The propagation of a wave in an initially uniformly distributed gas, is governed 
by the basic equations of mass, momentum, and thermal energy conservation 
conservation. With gas density p, at pressure p, and temperature T these equa- 
tion are: 

g + V.(pv)=0, (1) 
Dw 

p-^+Vf> = 0, (2) 

where v is the velocity of the wave in the gas, and Q(p,T) is the net heat gain 
per unit volume per unit time. It expresses the difference between the heating 
and cooling rates of the medium. Here A is the thermal conductivity coefficient, 
given by the ratio of the of specific heats 7 = c p /c v . 

For a gas at fixed volume and mass, the net heat gain is a function of temper- 
ature only. It therefore follows from equation (1) to (3) that instability occurs 
for the gaseous medium when ^ a h?r ) > (Field 1965). This condition is the 

Field stability criterion. In one spatial dimension the Field criterion implies 
that wave motion of heating or cooling interfaces must be highly subsonic. This 
means that we can neglect the velocity term in the momentum equation, hence 
obtain from equation (2): 

implying that pressure is constant within the medium. A change to Lagrange 
variable rj(x,t) defined by r/(x,t) = p(x' ,t)dx' , gives for the mass conser- 
vation equation (I): 

dt ' x= ~ pVl ^ 
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where v is the velocity in x-direction. This gives d/dx \ t = pd/drj \t&nd d/dt \ x = 
d/dt \rj —pud/drj \ t . The Lagrange derivative D/Dt can now be expressed as: 

at u +V T X '«= at ' (6) 

Using these results in the energy equation (3) we obtain, 

1 D P y pDf> 8 f^\ +Qi (7) 



7 — 1 Dt 7 — 1 p Dt dx \ dx 
which becomes 

7 pdp d ( dT\ 

In the subsonic case the equation of state is p = (TZ/ p) pT, it therefore follows 
that p is proportional to 1/T. Thus both A and Q are functions of temperature 

only. Let now the new time variable r = ^(7— l)p/7 (TZ/p) 2 ^j t, and denote 

C(T) = Q(T) (y—^J /T. Then equation (10) gives the governing equation for 
the evolution of the temperature front in a thermally conducting gas: 



dT_d_ ( A(T) 8T 



C(T). (9) 



3 Soliton Wave Analysis 

The conductivity function A(T) for a neutral ion gas is proportional to\/T. 
In the case of steadily moving wave front in 77 space, with velocity V, we set 
£ = c(i] — Kt). This gives in equation (12): 



, d 2 u r ^du 

c 2 u- 1 —+cK—+F{u)=Q, (10) 



where u = y/T, and F(u) = ^T -1 / 2 'C(T) is the net heating function. We now 
consider equation (13) for the case of linear, and quadratic net heating functions. 

If the net heating function is linear, we have from equation (10) the nonlinear 
wave equation: 

2 ld 2 u T ^du 

c 2 -—+cK—+Au + B = Q. 11 
u dt, 2 dt, 

We seek the solution to this equation such that is a function of the 

separable form, 

u'(0= a 2 -u(0 2 . (12) 

This implies that 

u"(0 = -2a 2 u(0 + 2u(0 3 . (13) 
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Therefore equation (11) becomes after substitution from equations (12) and 
(13), and collecting terms of the like powers in u: 

(2c 2 - cK )u 3 + Au 2 + (cKa 2 + B- 2a 2 c 2 )u = 0. (14) 

Equating the coefficients of u 2 , u 1 , and u° we get 

K = 2c, A = 0, B = 0. (15) 

This case correspond to the vanishing net heating function, and from equations 
(12) and (15) the solution can be written as 

u(0 =atanha(£ + £o), (16) 

where £o corresponds to the initial value of the travelling wave variable. Equa- 
tion (16) gives for the temperature 

T(r), t) = a 2 tanh 2 ac(r] - 2cr + £ /c). (17) 

Equation (17) implies that as £ — > ±oo, T — > T± QC = a 2 . The aysmptotic 
limit of temperature must correspond to the equilibrium state. We re-scale the 
temperature as T(r], r) = T± oc — T(r), t), so that equation (17) becomes 

T(v, t) = T ±00 (sech ac(f] - 2cr + ^/c)) 2 . (17a) 

This represents a free soliton wave solution when F(u) — 0. Physically this 
corresponds to the case when heating of the gas occurs at the same rate as 
cooling. The soliton amplitude in this case is a 2 and its speed is 2ac 2 . As 
£ — > ±oo, the soliton attains a fixed amplitude a 2 , whereas it retains its profile 
during propagation. Figure (1) shows plot of the soliton (17a) as function of 
variables 77 and r. 

For a quadratic heating function F(u) = Au 2 + Bu + C, we have by a similar 
procedure: 

(A + 2c 2 - cK)u 3 + Bu 2 + (cKa 2 + C- 2a 2 c 2 )u = 0. (18) 
Equating coefficients of the various powers of u, gives the determining equations, 
A + 2c 2 - cK = 0, cKa 2 + C- 2a 2 c 2 , B = 0, (19) 

we get 

c = ^ ± ^y/K 2 -8A, a = ±y/C/A, 5 = 0, (20) 

Thus corresponding to the net heating function Au 2 + C we have the soliton 

solution 

u(0 = ±y/C/Ataiih±y/C/A[c{ri - Kt) + ^0], (21) 

or in terms of the temperature variable 



T( V ,t) = - tanh 2 ±^ — [(K+ VK 2 -8A)( V - Kr) + 4&\. (22) 
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Re-scaling the temperature as T(r/, r) = T± OQ — T(r], r), where T± OQ = C/A, we 
obtain 



2 , . / C 



T(v,t)=T ±00 scch ± y — [(K + 2 - 8A) fa - if r) + 4£ ] . (22a) 

Equation (22a) shows that the soliton solution for a purely quadratic function 
has a strong amplitude dependence on the heating function parameters. In this 
way it represents a totally distinct soliton solution than that for the previous 
case. The thermal soliton (22a) has amplitude C/A, which increases as the 
parameter A decreases. Also the wave speed is proportional to J C/16A, which 
shows that the soliton wave moves faster if the quadratic term in the net heating 
function is adjusted to smaller values as compare to the constant terms in the 
heating function. Another interesting feature of the soliton solution (22a) is 
that the argument of the sech function involves the square root \JK 2 — 8 A, 
which can become imaginary when K 2 < 8A, thus exhibit drastically different 
behavior depending on the threshold value K = ±y/8A. 



4 Conclusions 

In this communication it is reported that thermally conducting gases in neutral 
ionized states contain thermal solitons that participate in the general process 
of energy transfer in the medium. For a given linear or quadratic net heating 
function these soliton waves can form as single solitary waves or, depending on 
the parameters of the heating function, can appear as multiple soliton waves. 
In the quadratic case, the net heating function parameters not only determine 
their amplitude but their speed in the gas as well. The existence of free solitons 
shows that there can be energy transfer via solitons even when the cooling of 
the gas occurs at the same rate as its heating. In this case the amplitude as 
well as the speed of the soliton wave depends on the equilibrium temperature 
of the gas. Since temperature is inversely propotional to the density of the gas, 
the soliton effects can be observed in density profile of a neutral ion gas, under 
appropriate conditions. 

Neutral ion gases are most prevalent in astrophysical and atmospheric con- 
ditions, such as the solar corona, interstellar medium, and galactic clusters, 
ionosphere, upper atmosphere of the Earth and planets. 

However the thermal solitons can also be observed under ordinary laboratory 
conditions. It is predicted here that for the vanishing and quadratic net heating 
functions, cases which are easier to be produced in a lab with adjustable heating 
parameters for neutral ion gases, the existence of these thermal solitons can be 
demonstrated. As far as the author knows this experimental proof yet remains 
to be given. 
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Figure Caption: 

Figure 1: A free soliton wave based on equation (17a), formed at equilibrium 
temperature T± OQ = 0.1, £o = 0, and c = 0.2. 
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This figure "fig 1. JPG" is available in "JPG" format from: 



http://arxiv.org/ps/1209.3077v3 



